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[i] Aerosol particles in the atmosphere are composed of multiple chemical species. The 
aerosol mixing state, which describes how chemical species are mixed at the single-particle 
level, provides critical information on microphysical characteristics that determine the 
interaction of aerosols with the climate system. The evaluation of mixing state has become the 
next challenge. This study uses aerosol time-of-flight mass spectrometry (ATOFMS) data and 
compares the results to those of the Goddard Institute for Space Studies modelE-MATRIX 
(Multiconfiguration Aerosol TRacker of mixing state) model, a global climate model that 
includes a detailed aerosol microphysical scheme. We use data from field campaigns that 
examine a variety of air mass regimens (urban, rural, and maritime). At all locations, polluted 
areas in California (Riverside, La Jolla, and Long Beach), a remote location in the Sierra 
Nevada Mountains (Sugar Pine) and observations from Jeju (South Korea), the majority of 
aerosol species are internally mixed. Coarse aerosol particles, those above 1 pm, are typically 
aged, such as coated dust or reacted sea-salt particles. Particles below 1 pm contain large 
fractions of organic material, internally mixed with sulfate and black carbon, and few external 
mixtures. We conclude that observations taken over multiple weeks characterize typical air 
mass types at a given location well; however, due to the instrumentation, we could not 
evaluate mass budgets. These results represent the first detailed comparison of single-particle 
mixing states in a global climate model with real-time single-particle mass spectrometry data, 
an important step in improving the representation of mixing state in global climate models. 
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1. Introduction 

[ 2 ] Atmospheric particles are rarely observed to be com- 
posed of a single chemical species [ Whiteaker et al., 2002], 
yet still the scientific community [IPCC, 2007] discusses 
the fate of sulfate aerosols and dust aerosols, among others, 
as though they exist independently. This terminology is 
most likely influenced by the way we have historically 
observed and modeled aerosol particles. The first mea- 
surements used filter-based techniques to record the mass 
of a chemical species [Malm et al., 2011], while most 
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models have simulated separate aerosol life cycles for 
each chemical species. 

[ 3 ] Details on the aerosol mixing state, such as the chemical 
components, water content, size, and shape (homogeneous 
mixed, core shell, etc.) of each particle, are of great importance 
for determining aerosol climate impacts. For example, a 
coated black carbon particle absorbs more radiation than its 
uncoated counterpart [Jacobson, 2006; Adachi et al., 2010], 
while a coated dust surface [Bauer et al., 2007] accelerates 
its wet removal and coagulated sulfate-organic carbon parti- 
cles provide number concentrations for cloud condensation 
nucleus (CCN) activation different from those counted 
separately [Bauer et al., 2010], Given the impacts that changes 
in mixing state can have on the climate properties of aerosol 
distribution, it is necessary to accurately model their properties 
with enough complexity to capture these processes. 

[ 4 ] To date, most measurements on the global scale are 
from space observations with remote sensing instruments 
(advanced very high resolution radiometer, Stratospheric 
Aerosol and Gas Experiment (SAGE), Total Ozone Mapping 
Spectrometer, Moderate Resolution Imaging Spectroradio- 
meter, Multiangle Imaging Spectroradiometer, Ozone Moni- 
toring Instrument, and Cloud-Aerosol Lidar and Infrared 
Pathfinder Satellite Observations (CALIOP)) on various satellite 
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platforms. This has led to a valuable source of global aero- 
sol maps that measure aerosol optical depth (AOD), absorp- 
tion AOD, and aerosol height, as well as additional aerosol 
quantities describing the total aerosol loading. Other long- 
term ground observations, such as the Interagency Moni- 
toring of Protected Visual Environments network in the 
United States and the European Monitoring and Evaluation 
Programme in the European Union, monitor time series 
of particle size and chemical mass loads. However, these 
remote and long-term measurements evaluate bulk aerosol 
distributions and, as the studies above have shown, 
knowledge of single-particle mixing state is necessary to 
evaluate aerosol impacts on climate. To investigate these 
processes, coordinated field campaigns provide in situ 
measurements over a shorter period of weeks to months 
and space with increasingly complex aerosol instrumenta- 
tion. Aerosol mass measurements, optics (Sun photometer, 
SP2), hygroscopicity, size distributions, differential mobil- 
ity particle sizer (DMPS), scanning mobility particle sizer 
(SMPS), and chemically resolved mass distributions (mass 
spectrometry such as AMS) have all provided insight into 
the complexity of aerosol particles on a bulk or average 
level. Single-particle measurements go a step further and 
allow for mixing state to be evaluated for individual parti- 
cles. In particular, online techniques such as particle analy- 
sis by laser mass spectrometry (PALMS) [ Murphy et al., 
1998], SPLAT [Zelenyuk et al., 2009], RSMS [Carson 
et aL, 1995], and aerosol time-of-flight mass spectrometry 
(ATOFMS) [Card et al., 1997] allow observation of the 
chemical composition of single particles with high statis- 
tics, which finally allows us to study the details of the 
mixing state of aerosol particles. 

[ 5 ] The first climate models that actively included aerosols 
contained the aerosol cycle by chemical species, beginning 
with the sulfur cycle [Schult et al., 1997] and then 
implementing additional chemical species (e.g., black carbon 
[Koch, 2001] and dust [Tegen and Miller, 1998]). However, 
the aerosol species were independent of one another. The 
inclusion of heterogeneous chemistry, particularly of the 
condensation of gaseous precursors on aerosol surfaces, in 
aerosol models [Dentener et al., 1996; Bauer, 2005] made 
aerosol transformation and removal dependent on other 
species and contained some consideration of mixing state. 
Aerosol microphysical models that resolve such processes as 
new particle formation, coagulation, condensation, and cloud 
activation do potentially consider many processes that lead 
to internally mixed aerosols. However, due to the need to 
keep aerosol simulations computationally efficient, com- 
promises had to be made when deciding on how to track this 
information. Most models are designed to keep track of the 
aerosol mass and number concentrations by using bin 
[Adams and Seinfeld, 2002; Spracklen et al., 2006] or modal 
[Whitby and McMurry, 1997; Easter et al., 2004] schemes. 
Even though model designs differ from one another, aerosols 
are commonly treated as either externally or internally 
homogenously mixed within one bin or mode. Exceptions 
are the sectional model by Jacobson [2002], which explicitly 
keeps track of aerosol mixing state, and the model used in this 
study, Goddard Institute for Space Studies (GISS) modelE- 
MATRIX (Multiconfiguration Aerosol TRacker of mixing 
state) [Bauer et al., 2008], a computationally efficient micro- 
physical aerosol model designed for climate applications. 


[6] The complexity of the problem and the difficulties in 
defining and observing aerosol mixing states make their 
evaluation an even greater challenge. This paper presents a 
first attempt at evaluating mixing states by comparing GISS 
modelE-MATRIX model findings to measurements made 
by an ATOFMS real-time single-particle mass spectrometer 
during five field campaigns. The campaigns took place in 
four locations in California and one location in South 
Korea during the years 2005-2009. 

2. Modeling Mixing State 

[ 7 ] The updated version of GISS modelE, which was 
prepared for GISS modelE contributions to the Coupled 
Model Intercomparison Project Phase 5 (CMIP5) archive, is 
used in this study. The model grid has horizontal and vertical 
resolutions (2° latitude x 2.5° longitude, 40 layers). In order 
to simulate similar tracer transport as observed during the 
field campaigns, model runs are nudged to the horizontal 
wind components of MERRA [Rienecker et al., 201 1], The 
model is run for the years 2000-2009 to cover the observa- 
tional periods of field campaign data used in this study. Sea 
surface temperatures are prescribed as boundary conditions 
and the model is fully coupled to the gas-phase chemistry 
scheme. Dust and seawater emissions are calculated interac- 
tively in the model, while anthropogenic and biomass 
burning emissions are from the CMIP5 emission inventory 
by Lamarque et al. [201 1]. 

[s] Aerosol processes in GISS modelE are calculated by 
the microphysical scheme MATRIX [Bauer et al., 2008, 
2010], MATRIX is based on the quadrature method of 
moments [Mcgraw, 1997], including two moments, number 
and mass, and one quadrature point. Carrying only two 
moments requires additional information about the shapes 
of the individual aerosol size distributions. We assume a 
lognormal distribution with constant width when calculating 
the initial size distributions, the conversion between aerosol 
mass and number concentration, emission distributions, 
coagulation rates, and aerosol optical properties. For each 
aerosol population (Table 1), defined by mixing state, the 
tracked variables are number concentration and mass concen- 
tration of sulfate, nitrate, ammonium, aerosol water, black 
carbon, organic carbon, mineral dust, and sea salt. 
MATRIX dynamics include new particle formation, particle 
emissions, gas-particle mass transfer, aerosol phase chemis- 
try, condensational growth, coagulation, and cloud activa- 
tion. An illustration of MATRIX size-resolved mixing state 
classes for distinctively different air mass types is presented 
in Appendix A. 

3. Observing Mixing State: Real-Time 
Single-Particle Mass Spectrometry 

[ 9 ] The aerodynamic size and chemical composition of 
individual particles between 0.1 and 3.0 pm were measured 
in real time using ATOFMS. The design and details of 
ATOFMS have been reported by Card et al. [1997], and a 
description is given in Appendix B. Single-particle size 
and mass spectral information were analyzed with YAADA 
1.2 (www.yaada.org), a data analysis toolkit for MATLAB 
6.5.1 (MathWorks, Inc.). Mass spectra were clustered using 
an adaptive-resonance-theory-based neural network algorithm 


9835 


BAUER ET AL.: AEROSOL MIXING STATE 


Table 1 . Definition of Mixing State in MATRIX: Population Definition, Geometric Mean Diameter at Emission Time Step, Name (as Used 
Hereafter), Chemical Composition, and Classification Compared to the ATOFMS Mixing State Classes 3 


Population Description 

Dg,E 

Symbol 

Constituents in Addition to S0 4 2 , N0 3 , NH 4 + , H 2 0 

ATOFMS Class 

Sulfate Aitken mode 

0.013 

AKK 


Sulfate/nitrate (other) 

Sulfate accumulation mode 

0.068 

ACC 


Sulfate/nitrate (other) 

Dust accumulation mode (<5% inorganic) 

0.58 

DD1 

Dust 

Fresh dust 

Dust accumulation mode (>5% inorganic) 


DS1 

Dust 

Aged dust 

Dust coarse mode (<5% inorganic) 

5.4 

DD2 

Dust 

Fresh dust 

Dust coarse mode (>5% inorganic) 


DS2 

Dust 

Aged dust 

Sea salt accumulation mode 

0.37 

SSA 

Sea salt 

Salts 

Sea salt coarse mode 

3.93 

SSC 

Sea salt 

Salts 

Organic carbon (OC) 

0.03 

OCC 

OC 

OC 

Black carbon (BC) (<5% inorganic) 

0.03 

BC1 

BC 

Fresh EC 

BC (5%-20% inorganic) 


BC2 

BC 

Aged EC 

BC (>20% inorganic) 


BC3 

BC 

Aged EC 

BC-OC 

0.037 

BOC 

BC, OC 

EC-OC 

BC-sulfate (formed by coagulation) 


BCS 

BC 

EC-sulfate 

BC -mineral dust 


DBC 

Dust, BC 

Dust-EC 

Mixed 


MXX 

Dust, sea salt, BC, OC 

SS-EC, dust aged, salt, others 


a For comparison of the available data, the modeled black carbon is compared to the observed elemental carbon. EC and BC are used synonymously. 


(ART-2a) at a vigilance factor of 0.8 and regrouped at a vigi- 
lance factor of 0.85 [Song et al . , 1999]. General particle types 
are defined by the characteristic chemical species or possible 
source to simplify the naming scheme. ATOFMS counts were 
scaled to number concentrations using aerodynamic particle 
sizer (APS) size-resolved number concentrations, a method 
shown previously to yield quantitative mass concentrations 
[Qin et al., 2006]. For the Gosan study, SMPS data were used 
to expand the lower limit of the scaling by applying number 
fractions to the SMPS size distributions after correcting for 
density, an imperfect method necessary for comparisons near 
the mode of the number size distribution. 


4. Comparison Between ATOFMS and MATRIX 

[ 10 ] The following comparisons demonstrate the ability of 
MATRIX to reproduce the size-resolved chemical composi- 
tion of aerosols in a number of environments: coastal marine, 
aged urban, fresh urban, and remote continental. 

[n] ATOFMS data are used for a number of field 
campaigns between 2005 and 2009, as listed in Table 2. 
For this comparison, we have sampled the observations over 
the entire duration of each field campaign, which lasted 
between 9 and 45 days (average of 22 days). The model 
was sampled for the same time episodes. To allow compari- 
sons between ATOFMS and MATRIX, the data sets had to 
be merged into similar mixing state definitions. The grouping 


of MATRIX populations into the ATOFMS classes is shown 
in Table 1. The general ATOFMS classes of fresh elemental 
carbon (EC) (soot), elemental carbon/organic carbon 
(ECOC), organic carbon (OC), sea salt, mineral dust, and 
mixtures account for >85% of particles from 0.1 to 3.0 pm 
in studies across the Northern Hemisphere. To bring the 
model categories and ATOFMS classes into agreement, 
ATOFMS clusters were recategorized based on the following 
distinctions: Fresh EC went into BC1, ECOC went into EC- 
OC, biomass burning and OC went into OCC, sea salt went 
into SSA or SSC (depending on size), fresh dust with 
minimal markers for aging (i.e., nitrate) went into DD1 and 
DD2 as a function of size, and aged dust with nitrate went 
into DS1 and DS2 as a function of size. Spectra with no 
positive signal and intense sulfate peaks went into ACC, 
mixtures of EC with dust into DBC, and mixtures of carbona- 
ceous and mechanically generated particles (i.e., EC with SS/ 
dust) into MXX, but none of these categories were significant 
number fractions at any size sampled. Modeled particle sizes 
have been converted into aerodynamic size by taking into 
account the density of the aerosol mixture. 

4.1. La Jolla, California 

[ 12 ] Measurements made at the end of the Scripps 
Institution of Oceanography (SIO) pier (~300 m offshore) 
in 2006 provide a representative example of a coastal marine 
location with urban influences, which due to the diurnal wind 


Table 2. ATOFMS Study Details 


Study Location 

Abbreviation 

Study 

Latitude 

Longitude 

Period 

References 

Riverside, California, 
USA 

SOAR LABMS1 

Study of Organic Aerosols 
in Riverside Los Angeles 
Mobile Study 1 

33.97 

-117.32 

30Jul. to 15 
Aug. 2005, 29 
Aug. to 7 Sep. 2007 

Moffet et al. [2008], 
Ault et al. [2009], 
Gaston et al. [2010] 

La Jolla (San Diego), California, 
USA 

SIO Pier 

Scripps Oceanography 
Pier Study 

32.87 

-117.26 

18 Aug. 
to 1 Oct. 2006 

Ault et al. [2009] 

Gosan, Jeju, South Korea 

PACDEX 

Pacific Dust Experiment 

33.29 

126.18 

13 Apr. 

to 15 May 2007 

Ault and Pomeroy [2012] 

Los Angeles, California, USA 

LABMS2 

Los Angeles Mobile Study 2 

33.74 

-118.27 

17-26 Nov. 2007 

Ault et al. [2010] 

Sugar Pine, California, USA 

Calwater 

Calwater Early Start 

39.13 

-120.80 

21 Feb. 

to 1 1 Mar. 2009 

Ault et al. [2011], 
Creamean et al. [201 1] 
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patterns did not experience significant stagnation [Ault et al., 
2009], Figure la shows size-resolved number concentration 
data averaged over the course of the study. The supermicron 
size range is dominated by sea-salt particles, with the mode 
peaking at roughly 1.1 pm. Below 1 pm, a transition is 
observed toward predominantly carbonaceous aerosols, 
which dominate the overall number concentration below 
0.6 pm. The diameter of the transition from >50% mechani- 
cally generated aerosols (i.e., sea salt and mineral dust) to 
>50% carbonaceous aerosols (soot, organic carbon, etc.), 
ri carb , has been identified as a useful parameter for evaluating 
how important sources other than sea salt and dust are on a 
site. For the SIO study, this value occurs at 0.64 pm. In this 
environment, the vast majority of carbonaceous particles 
contain elemental carbon (or soot) mixed with such second- 
ary species as sulfate, nitrate, organic carbon, and water 
and there is no significant externally mixed population of 
organic carbon particles above 0.5 pm. 

[b] Comparison between ATOFMS and MATRIX (Figure 1) 
shows a remarkable agreement between the model and the 
observations, in terms of both number concentration and 
single-particle composition. Among the key features in 
agreement between the model and the measurements are 
that the supermicron sea-salt mode peaks at 1 . 1 pm and that 
carbonaceous species dominate below 0.6 pm. There are 
also a few noteworthy differences that have been observed 
between the model and the measurements, including a lack 
of mineral dust in ATOFMS measurements compared to 
MATRIX predictions, the relative contribution of purely 
organic carbon versus aged soot to the carbonaceous mode 
below 0.6 pm, and a lack of fresh EC. The modeled dust 
originated from Mexico and was transported to this location 
carrying sulfate, ammonia, and nitrate coatings on its sur- 
face. This discrepancy may be due in part to the comparison 
of a grid cell size area having significant impacts from dust 
to a site that has minimal impact from dust except under 
strong offshore wind conditions (Santa Ana events) 
[Guazzotti et al., 2001]. Mixing state classes of particles 
sized between 0.5 and 1 pm are dominated by mixtures 
incorporating organics. The largest fraction is ECOC, 
followed by OC and partial contributions of EC-sulfate, 
which is larger in the observations than in MATRIX. The 
lack of fresh EC (low sulfate, nitrate, or organic carbon 
loadings) is simulated at all locations shown in Figure 1 
and likely due to local sources not captured by the model 
or averaged out over the size of the grid cell. Capturing fresh 
EC and understanding its aging in models and measure- 
ments at the same sites are an ongoing challenge that calls 
for further efforts between modelers and experimentalists. 
We will revisit this topic later for the Jeju case, where we 
can compare mixing state information between ATOFMS 
and MATRIX for smaller particle sizes. 

4.2. Riverside, California 

[ 14 ] In contrast to San Diego, where diurnal winds serve to 
prevent stagnation, Riverside, California, is one of the most 
prominent examples of a location where air masses stagnate, 
allowing for extensive formation of secondary organic aero- 
sols [Docherty et al., 2008; Pratt and Prather, 2009; Qin 
etal., 2012]. Data from the summer of 2005 (Figure lb) show 
a substantially different single-particle composition, with the 
organic carbon particles representing the majority of particles 


below 1 .4 pm, with the rest coming from fresh EC and ECOC 
[Qin et al., 2012], The d cm b transition to the majority being 
mechanically generated particles occurs at 1 .4 pm, a much 
higher value than that for San Diego, indicating that 
supermicron carbonaceous particles are a significant fraction 
of supermicron mass in stagnant locations. In contrast to 
2005, the data from Riverside in 2007 (Figure lc) show a 
period with much lower PM concentrations and a greatly 
reduced contribution of carbonaceous particles above 1 pm 
[Ault et al., 2010]. A sea-salt mode that peaks at 1.2 pm, 
observed in both San Diego and Riverside during 2005, 
dominates the supermicron micron number fraction. Due to 
the high winds that resulted in less stagnant conditions, 
mineral dust is 5 %— 1 5 % of the number fraction between 
0.5 and 3.0 pm. Below 1 pm, internally mixed ECOC is the 
greatest contributor to the number concentration, with 
smaller contributions from fresh EC and OC particles. 
These studies show that while the intensity of the carbona- 
ceous particle dominating the number size distribution varies 
in Riverside, the mechanically generated sea-salt and dust 
concentrations are reasonably consistent under different 
meteorological conditions. 

[ 15 ] Very similar distributions from ATOFMS and 
MATRIX are shown for Riverside (Figures lb and lc). 
During the 2005 field observational period, MATRIX shows 
a smaller number concentration, but particles below 1 pm are 
primarily ECOC particles, OC, and a small fraction of EC- 
sulfate, with good agreement as to their relative contribu- 
tions. Similarly, the model captures a transition to sea salt 
and dust above 1 pm, although the diameter at which this 
occurs is smaller than that for the ATOFMS measurements. 
Additionally, MATRIX finds a significant contribution of 
sea salt and mineral dust down to 0.5 pm that is not observed 
in the ATOFMS measurements. During the 2007 period, 
MATRIX captured smaller contributions of pollution 
aerosols and modes of sea salt and mineral dust. MATRIX 
does show a slightly higher fraction of dust, while the 
partitioning between fresh dust and aged dust shows more 
fresh dust than the measurements indicate. Interestingly, 
both examples show a small fraction of sulfate aerosols 
(marked in yellow). Overall, MATRIX does well at 
representing the single-particle composition of the number 
size distribution for a location that has conditions ranging 
from extremely stagnate and polluted to relatively clean un- 
der high wind conditions. 

4.3. Long Beach, California 

[ 16 ] A challenging location to sync measurements and 
models is a region with multiple intense and rapidly 
evolving sources, such as the region surrounding the ports 
of Los Angeles and Long Beach [Krudysz et al., 2009]. 
Measurements on Terminal Island in the fall of 2007 repre- 
sent an extreme case, with diurnal wind patterns preventing 
stagnation but extensive sources from shipping, trucking, 
and industry (i.e., oil refineries) leading to extremely 
polluted conditions [Minguillon et al., 2008; Krudysz 
et al., 2009; Ault et al., 2010], Here ATOFMS data show 
high concentrations of fresh EC and ECOC dominating 
the number size distribution below 1 .4 pm and that, unlike 
in Riverside, purely organic carbon particles make up less 
than 5% of aerosols in all bins above 0.5 pm. The sea-salt 
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a) 


La Jolla, San Diego, September 2006 



b) 


Aerodynamic diameter [pm] 

Riverside, CA July-August 2005 



1.0 1.5 2.0 

Aerodynamic diameter [urn] 




1.0 1.5 2.0 2.5 

Aerodynamic diameter [urn] 

Sugar Pine, Lake Tahoe, Feb/March 2009 
l a.oot 



Aerodynamic diameter [urn] 


MATRIX 



MATRIX 



MATRIX 



Aerodynamic diameter [pm] 

MATRIX 



MATRIX 



Aerodynamic diameter [pm] 


EC fresh 

EC-OC 

EC-Sulfate 

OC 

Salts 

DustSupFresh 

DustSupAged 

Dust-EC 

SS-EC 

Other 


Figure 1. Size-resolved number concentrations of color-coded mixing state classes are shown for particle 
size bins between 0.5 and 3 pm. The inserted graphs show the same distributions but on a y-log axis. 
ATOFMS (left) and MATRIX (right column) results are shown for five field campaigns: (a) one in La 
Jolla, (b and c) two in Riverside, and one each in (d) Long Beach and (e) Sugar Pine. 


mode observed in the previous locations is observed again 
at -1.3 pm, but — as with Riverside in 2005 — the anthro- 
pogenic particles dominate the distribution. 

[n] The source site represented is very spatially homoge- 
neous. The model did not capture these features and 
involved more general background conditions, similar to 
those observed for Riverside in 2005. Given the large num- 
ber of high-intensity sources, the ports of Los Angeles and 


Long Beach represent a site influenced by local pollution 
that cannot be captured by the general circulation model 
grid boxes. 

4.4. Sugar Pine, California 

[is] A second challenging location with extremely low par- 
ticle concentrations is a remote continental site. Measurements 
at the Sugar Pine Reservoir in Forest Hill, California, a remote 
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105 110 115 120 125 130 135 


Longitude (Degrees East) 

Figure 2. Size-resolved number concentrations of mixing 
state classes for particles between 0.07 and 0.5 pm in Jeju, 
South Korea, for ATOFMS (left) and MATRIX (right). The 
inserted boxes show the entire range of measured size distri- 
bution on a log scale. Averaged mixing state conditions are 
shown for four air mass types originating from (a) a marine 
site, (b) Eastern China, (c) Korea, and (d) the Gobi desert. 

rural site in the Sierra Nevada Mountains (elevation, -1064 m 
ASL) [Ault et al., 2011; Creamean et al., 2011], provide an 
ideal test case for evaluating MATRIX in such a location. 
In comparison to Long Beach, which routinely has number 
concentrations >10,000 #/cm 3 [Krudysz et al., 2009], num- 
ber concentrations at Sugar Pine averaged 240 #/cm 3 
[Creamean et al., 2011], roughly 2 orders of magnitude 
lower. The number size distributions at Sugar Pine were 
frequently dominated by new particle formation events, 
meaning that nearly all of the number concentrations were 
below 0.15 pm [Creamean et al., 2011]. Thus, the very 
few particles >0.5 pm led to the low concentrations shown 
in Figure le. In the supermicron range, much of the particle 
concentration is composed of sea salt and a mix of fresh dust 
and aged dust, while below 1 pm, ECOC, OC, and EC- 
sulfate comprise the majority of the particles. The sea-salt 
mode results from long-range transport to the site [Ault 
et al., 2011] and is at the same diameter as observed in the 
other locations. A shift is observed from equal parts aged 
dust and fresh dust in the supermicron to primarily fresh 
dust in the submicron. 

[ 19 ] For the remote continental site, MATRIX does a good 
overall job of capturing the different particle compositions 
and their sizes. In the supermicron mode, the composition 
is dominated by dust and sea salt in similar proportions to 
the ATOFMS measurements, although all dust is observed 
as aged due to the model only considering desert dust 
sources. Below 1 pm, contributions from ECOC, EC-sulfate, 
and OC are shown, with little fresh EC matching the 
ATOFMS measurements, at roughly similar relative propor- 
tions. The shape of the dust distribution in MATRIX matches 
the ATOFMS observations. The challenge of grid size is 


observed again, although in this case with much higher 
concentrations in MATRIX than in measurements versus 
the opposite trend observed at the ports of Los Angeles and 
Long Beach. This is logical as other sources exist in the grid 
cell (such as 1-80 and several small cities) and Sugar Pine was 
chosen to be as pristine a site as possible. However, from a 
fractional perspective, MATRIX does well for a very 
challenging site. 

4.5. Jeju, South Korea 

[ 20 ] To test the agreement between MATRIX and single- 
particle measurements closer to the peak of the number size 
distribution, data from Gosan, South Korea, in 2007 are 
shown in Figure 2. The island is impacted by air masses 
containing particles from a wide range of source regions 
(Figure 3), including heavily polluted air masses from 
Eastern China, less polluted air masses from the mainland 
of Korea, air masses from the Gobi desert, and marine air 
masses. The number size distributions peaked between 
0.067 and 0.111 pm depending on the air mass types. 
ATOFMS data down to 0.1 pm give a good indication as 
to the different single-particle compositions likely to be 
present at the mode of the distribution. It is interesting to 
observe that mixing state classes do not dramatically change 
between the four air mass cases: marine, East China, Korea, 
and Gobi desert. Typically, the composition includes a 
small fraction of dust (except in the marine case), a large 
contribution from internally mixed OC, having EC and 
sulfate, and a substantial contribution of fresh EC (due to 
local sources). The <r/ carb is observed at roughly 1 pm, above 
which sea salt and dust become dominant. Of particular note 
is that the greatest difference between the air mass types is 
seen in the absolute number concentrations, which are 
highest when the air mass has spent time over Eastern 
China and lowest for the marine back trajectory. MATRIX 
does well at capturing the differences in number concentra- 
tion between the four cases. However, even with these wide 
ranges in number concentration, the composition of 
particles <0.5 pm does not change significantly, which is 
important for how marine regions are modeled generally. 
This is caused by the fact that Gosan experiences relatively 
aged air masses. It is important to note the lack of pure 
OC particles, which are often assumed to be large contribu- 
tors to the submicron aerosol in remote locations. ATOFMS 
and MATRIX show a significant contribution of fresh 
EC. Fresh EC, which has been detected by ATOFMS 
(Figure 1), is now finally visible in MATRIX as well 
because we are comparing smaller aerosol size classes in 
this case while in a region that is generally very high in 
black carbon emissions. Fresh EC coagulates and attracts 
coatings quickly, allowing fresh EC in the model only to 
be visible around its emission sizes. Larger EC that has 
already aged will always appear in mixed EC size classes. 
Thus, the apparent lack of consistency between ATOFMS 
and MATRIX measurements for the studies in Figure 1 is 
likely definitional and does not represent significantly 
different results. 

[ 21 ] Figure 4 shows the total number concentrations for 
particles between 0.5 and 3 pm for all six field campaigns. 
For number concentrations above 5 cm -3 , the model and 
observations show similar scatter. For concentrations below 
5 cm -3 , mostly representing supermicron aerosols, the 
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Figure 3. Back trajectories of air mass types as observed in Jeju, South Korea, between 12 April and 15 
May 2007. 


model tends to underestimate EC-sea salt and EC-dust 
mixtures, while aged dust, fresh dust, and salts in this 
concentration range are overestimated. Dust and sea salt 
are the aerosol types that spread over all size ranges. 
MATRIX has to resolve the size distributions and mixing 
state of dust with six populations and those of sea salt with 
three. These resolutions might not be enough to represent all 
possible size and mixing state options precisely. Furthermore, 
we only emit dust and sea salt into two size classes, which is a 
simplified assumption. Lastly, the definitions of mixing state 


classes between the measurements and the model are not iden- 
tical and are only approximated. 

5. Discussion 

[ 22 ] The microphysical details that are necessary to simu- 
late realistic aerosol mixing states require the model to have 
good skills in general aerosol simulation, such as aerosol 
mass by chemical species, number concentrations, and 
aerosol size. From prior [Bauer et al., 2008, 2010, 2013] 
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Figure 4. Scatter plot between model and ATOFMS total 
number concentrations (summed over 0. 5-3.0 pm) per 
mixing state in cm -3 . All stations from Figure 1 and Jeju 
are presented. The color coding indicates the mixing state 
classes. 

and ongoing evaluations through the AeroCom initiative 
(http://aerocom.met.no), we have learned the strengths and 
weaknesses of MATRIX relative to those of other global 
aerosol models. Therefore, we are pleased to leam from this 
study that comparisons between a coarse global climate 
model and detailed field measurements at single locations 
are possible. Averaging the field campaign data over the 
duration of the campaign, typically a couple of weeks, 
allowed for the characterization of typical air masses and 
mixing state classes at each location. Further nudging the 
model with reanalysis data allowed the comparison of similar 
air mass types that have been transported to the observed 
location. Although not discussed in this paper, conducting a 
similar comparison using a climatological model run (10 year 
mean conditions of a free-running, nonnudged model exper- 
iment) found mostly similar results as presented in this study. 
Hence, we conclude that observations taken over a period of 
1 month sufficiently characterize typical air mass types at a 
given location. 

[ 23 ] This comparison has demonstrated that at the four 
locations in California, three in polluted areas (Riverside, 
La Jolla, and Long Beach) and one in a remote location in 
Sugar Pine, the great majority of aerosol species are inter- 
nally mixed. Coarse aerosol particles, above 1 pm, tend to 
be mechanically produced and aged, such as coated mineral 
dust or heterogeneously reacted sea-salt particles. Particles 
below 1 pm contain large fractions of organic material and 
appear in mixing states that include organics such as EC- 
OC and OC. No externally mixed sulfate or nitrate was 
observed or simulated in the presented examples in particles 
larger than 0.5 pm. Observations in Jeju, South Korea, 
allowed us to compare particles down to 0.1 pm. At these 
smaller sizes, MATRIX could detect more fresh EC, in 
agreement with the measurements. Overall, the ATOFMS 
definition of fresh EC (carbon cluster peaks in the positive 
and negative mass spectra, with greater intensity than OC 
or sulfate peaks) and MATRIX definition of fresh EC 



(<5% inorganic and no OC) were not in complete agreement 
(as seen in Figure 3) but converged at lower sizes where more 
fresh EC is present in the model (Figure 2). Finding defini- 
tions that can help lead to convergence is both challenging 
and necessary given the importance of representing EC/soot 
mixing state properly for climate impacts [ Jacobson , 2002; 
Riemer et al., 2009]. 

[ 24 ] This study is limited by the fact that we only compared 
number concentrations of mixing state classes above 0.5 pm 
(and in the Korea case, above 0.07 pm) diameters. Particles 
above 0.5 pm have already aged significantly through 
condensation and coagulation processes. Also, we did not 
compare mass budgets in this study. As a logical next step, 
we would like to evaluate the mixing states of fine and coarse 
aerosol concentrations to close mass and number budgets. A 
multi-instrument analysis is most suitable to deal with a 
complex aerosol property such as mixing state. For example, 
the combination of state-of-the-art single-particle mass spec- 
trometer instruments such as a single-particle mass spectrom- 
eter (SPLAT II [Zelenyuk et al., 2009]), the particle analysis 
by laser mass spectrometer (PALMS [. Murphy et al., 1998]), 
a single-particle soot photometer (SP2 [Schwarz et al., 
2008]), and the aerosol time-of-flight mass spectrometer 
(ATOFMS [Pratt et al., 2009]) used in this study can provide 
a more complete picture of the different particle types and 
mixing states present in a population of aerosols. Aerosol 
properties such as density, optical properties, shape, number 
concentrations, and size distributions [Zaveri et al., 2012] as 
well as data from other sources to determine composition as a 
function of hygroscopicity and CCN activity, among others, 
will further provide insight into the aerosol mixing state and 
its relevance for climate simulations. 

[ 25 ] Resolving mixing state in climate models, especially 
when attempting to answer questions regarding which emis- 
sion source is responsible for air pollution or climate impacts, 
is of great importance. Instruments such as ATOFMS pro- 
vide unique details for identifying the potential sources of 
aerosols (e.g., vehicular emissions and biomass burning) as 
well as their degrees of aging. The combination of observa- 
tions and models that can resolve such chemical details is 
required for advancing climate and policy studies. 

Appendix A: Examples of Mixing 
State Representation 

[ 26 ] In order to illustrate the definitions of MATRIX 
mixing state classes, Figure Al presents size-resolved aero- 
sol population examples of different air mass types. Four ex- 
amples of annual mean mixing state conditions at the surface 
level are presented for marine (close to Hawaii), desert 
(Sahara, Mali), biomass burning (Brazil), and urban pollution 
(Beijing) grid cells. The size-resolved fraction of the chemi- 
cal composition is given for the seven components of the 
aerosol populations. Aerosol water, which is not shown in 
this plot, contributes to 79%, 1%, 31%, and 77% of the total 
aerosol mass in the chosen examples, respectively. Small 
particles, below 0.1 pm, in the remote marine location are 
composed of sulfate and ammonium, represented in the 
Aitken and accumulation (AKK and ACC) populations. Sea 
salt is the dominant chemical species of particles larger than 
0.1 pm and represented in the SSA, SSC, and MXX 
populations. MXX at this location mostly includes sea salt 
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Figure Al. Aerosol mixing state samples of various air mass types: remote marine (Hawaii), desert 
(Mali), biomass burning (Gabon, Brazil), and urban (Beijing). Size distributions (geometric diameter in 
micrometers) of mixing state fractions are given per chemical species (sulfate, nitrate, ammonium, BC, 
OC dust, and sea salt) and aerosol population (population acronyms are defined in Table 1). 


and a small fraction of dust. Small particles in the desert 
example are composed of sulfate, ammonia, nitrate, and 
organic carbon, with a small fraction of black carbon. 
Larger particles are completely dominated by dust minerals. 


The larger particles are dominated by different mixing states 
of desert dust. These include small and large dust classes with 
thin (DD1 and DD2) and thick (DS1 and DS2) coatings, a 
small fraction of BC-dust mixtures, and, finally, MXX, 
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which in this location primarily contains dust but with larger 
coatings than those that can be accommodated in DS1 or 
DS2. The location picked in Brazil experiences frequent bio- 
mass burning events. This leads to a high concentration of 
OC in the aerosol mix resulting in OCC, BC1, and BOC 
populations for aerosols of sizes between 0.1 and 1 pm. 
Coarse particles are composed of a mix of dust and sea salt 
with small fractions of inorganic species. Beijing accommo- 
dates larger fractions of carbonaceous aerosols compared to 
the previous examples and the larger particles are dominated 
by desert dust. The higher BC and OC concentrations lead to 
the spread of OCC, BOC, and BCS into the small aerosol size 
classes (0.1 pm). Additionally, nearly pure BC, population 
BC 1 , can be observed. 

[ 27 ] These examples are given to visualize how mixing 
state is represented in the model, over the size classes we 
resolve, from 0.03 to 20 pm. The mixing state comparisons 
between the measurements and the model presented in this 
paper only allow comparing the size classes that can be 
observed by ATOFMS, which in most cases are coarse aero- 
sol fractions above 0.5 pm and in one case cover particles 
from 0.07 to 3 pm. 

Appendix B: Aerosol Time-of-Flight 
Mass Spectrometry 

[ 28 ] Particles are introduced into the ATOFMS through a 
converging nozzle into a differentially pumped vacuum cham- 
ber where they are accelerated to a terminal velocity [Card 
et al., 1997], The particles then pass through two continuous 
wave lasers (diode-pumped Nd:YAG lasers operating at 
532 nm) located 6 cm apart. The aerodynamic diameter of 
each particle is determined through calibration with polysty- 
rene latex spheres of known aerodynamic size, shape, and 
density. The sized particles are then desorbed and ionized by 
irradiation from a Q-switched Nd : YAG laser (266 nm, 1.2— 
1 .4 mJ/pulse) that is triggered when the particle is in the center 
of the mass spectrometer ion source region. Positive and neg- 
ative ions from each particle are detected using a dual-polarity 
reflectron time-of-flight mass spectrometer. These labels do 
not reflect all of the species present within a specific particle 
type but do reflect the most intense ion peaks. 

[ 29 ] The scaling to number concentrations is done as fol- 
lows: ATOFMS counts are binned into sizes that correspond 
to the aerodynamic particle sizer (APS) bins (0.5-2. 5 pm) 
and scaling factors for each hour and size bin are determined 
to account for ATOFMS transmission biases and busy time. 
Each bin is then scaled to give number concentrations. 
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